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Abstract 



We show that the equilibrium positions of the Ruijsenaars- Schneider- van Diejen 
systems with the trigonometric potential are given by the zeros of the Askey- Wilson 
polynomials with five parameters. The corresponding single particle quantum version, 
which is a typical example of "discrete" quantum mechanical systems with a (/-shift 
type kinetic term, is shape invariant and the eigenfunctions are the Askey-Wilson 
polynomials. This is an extension of our previous study which established the 

"discrete analogue" of the well-known fact; The equilibrium positions of the Calogero 
systems are described by the Hermite and Laguerre polynomials, whereas the corre- 
sponding single particle quantum versions are shape invariant and the eigenfunctions 
are the Hermite and Laguerre polynomials. 



1 Introduction 



The Calogero- Sutherland systems [3] and their integrable deformation called the Ruijsenaars- 
Schneider-van Diejen systems jU |H] have many attractive features at both classical and 
quantum mechanical levels. In our recent paper Q], the equilibrium positions of the 
classical Ruijsenaars-Schneider-van Diejen systems were studied. The equilibrium positions 
of the Calogero-Sutherland systems are described by the zeros of the classical orthogonal 
polynomials, the Hermite, Laguerre, Chebyshev, Legendre, Gegenbauer and Jacobi polyno- 
mials pJIBllII]- Since the Ruijsenaars-Schneider-van Diejen systems are deformation of the 
Calogero-Sutherland systems, it is expected that the equilibrium positions of the Ruijsenaars- 
Schneider-van Diejen systems are described by some deformation of these classical orthogonal 
polynomials. This is indeed the case and we obtained the deformed Hermite, Laguerre and 
Jacobi polynomials pQ. These deformed orthogonal polynomials fit in the Askey-scheme of 
the hypergeometric orthogonal polynomials [TTH lllj: (i) rational potential cases: One and 
two parameter deformation of the Hermite polynomials are a special case of the Meixner- 
Pollaczek polynomial and a special case of the continuous Hahn polynomial, and two and 
three parameter deformation of the Laguerre polynomials are the continuous dual Hahn poly- 
nomial and the Wilson polynomial, (ii) trigonometric potential cases: Several one parameter 
deformation of the Jacobi polynomials are special cases of the Askey- Wilson polynomial. 
The Askey- Wilson polynomial has five parameters |T2], but the deformed Jacobi polynomi- 
als obtained in |]Q have only three parameters. A natural question arises; Find (integrable) 
multi-particle systems whose equilibrium positions are described by the Askey- Wilson poly- 
nomials with five parameters. 

Shape invariance is an important ingredient of many exactly solvable quantum mechanics 
[121 CH HH] • In another recent paper of ours [2], the shape invariance of "discrete" quantum 
mechanical single particle systems, whose kinetic term causes a shift of the coordinate in 
the imaginary direction, are discussed. The eigenfunctions of these shape invariant systems 
are a special case of the Meixner-Pollaczek polynomial, a special case of the continuous 
Hahn polynomial, the continuous dual Hahn polynomial and the Wilson polynomial. These 
polynomials have all appeared in the above discussion about the equilibrium positions, in 
which we have one more polynomial, the Askey- Wilson polynomial. This gives the second 
question; Are the quantum mechanical single particle systems, whose eigenstates are the 
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Askey- Wilson polynomial, shape invariant or not ? 

We will answer the above two questions in this paper. The answer to the first question 
is found by the same method given in pQ, i.e. numerical analysis, functional equation and 
three-term recurrence. The second question is answered affirmatively by using the properties 
of the Askey- Wilson polynomials and similar discussion in [S] with replacement of a shift 
operator by a g-shift operator. 

This article is organized as follows. In section 2, the classical equilibria of the Ruijsenaars- 
Schneider-van Diejen system are studied. For a suitable choice of the elementary potential 
functions, the equilibrium positions are given by the zeros of the Askey- Wilson polynomials 
with five parameters. In section 3 we discuss the shape invariance of "discrete" quantum 
mechanical single particle systems with a g-shift type kinetic term. After discussing general 
theory, we present an explicit example of such shape invariant systems, in which eigenstates 
are described by the Askey- Wilson polynomials. The final section is for a summary and 
comments. 

2 Multi-particle Systems: Equilibrium Positions 

Let us consider the equilibrium positions of the classical Ruijsenaars-Schneider-van Diejen 
systems with the trigonometric potential [UE], which are integrable deformation of the cel- 
ebrated Calogero- Sutherland systems of exactly solvable multi-particle quantum mechanics. 
Its classical Hamiltonian corresponding to the BC root system is the following [Hj: 

H( P , g)=E( cosh ^' yJ v MVjb) ~ « fad + > W 

3=1 ^ ' 

n 

Vj(q) = w(qj) Yl v(qj - q k ) v(qj + q k ) (j = 1, . . . , n), (2) 

fe=i 

Mz-ino) (3) 

smx 

W M - sin ( x ~ igi ^ sin ^ ~ cos ( x ~ ig ^> cos ( x ~ (4) 
sinx sin a; cosx cosx 

where q = t (qi, ■ ■ ■ ,q n ) and p = t (pi, ■ ■ ■ ,p n ) are the coordinates and conjugate momenta, 
and Qj (j = 0, . . . ,4) are the real positive coupling constants. The potentials Vj and V* are 
complex conjugate of each other. Our convention of a complex conjugate function is the 
following: for an arbitrary function f(x) = J2 n a n xn ( a n & C), we define f*(x) = J2 n a n xn - 
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Here c* is the complex conjugation of a number c G C. Note that f*(x) is not the complex 
conjugation of f(x), (f(x))* = f*(x*). This is relevant for considering complex variables in 
section 01 The equilibrium positions p = 0, q = q are determined by the condition 

V j (q) = V*(q)>0 (j = l,2,...,n). (5) 

This equation without inequality is rewritten in the Bethe ansatz like equation 

tt v(qj - Qk)v(qj + qk) _ w*(%) 

f£v*%-q k )v*(q j + q k ) w{qj) 

Note that q\j — 0, ~ is excluded in (JSJ) but allowed in (JHJ). 

By the same method given in {Tj (numerical analysis, functional equation and three- 
term recurrence), we can show that the equilibrium positions {q~j} are given by the zeros of 
the Askey- Wilson polynomial (we follow the notation of Koekoek and Swarttouw [TT] ) 
p n (x; a, b, c, d\q) oc nj=i( x — cos2g,) with the following parameters, 1 

q = e - 290 , (a,b jCj d) = (e- 29l ,e- 292 ,-e- 293 ,-e- 294 ). (7) 

The outline of derivation is as follows. The functional equation for f(x) = YYj=i ( x ~ cos ^Qj) 
is the same as eq.(A.4)(e = —1) in 1\ with 2 

h(x) = (*±=*x - i(l + 2li^)v / T^ 2 ") 

x (^(1 + x) + z(l + sli±2)v / T^ 2 ") (^(1 + x) + i(l + K^)Vl^) 

x (^(1 -x)- i(l - ^)VT^ 2 ) (^(1 -x)- z(l - ^)v / I^), 

(8) 

(1 + a)(l + b)(l - c)(l - d)(l + g)(l - x 2 ) 
8abcdq n+1 

x ^4g(l + abcdq 2n ~ x )x 2 - 2q n (abc + abd + acd + bed + (a + b + c + d)q)x 
- (1 + q) (l + q - q n+1 - (ab + ac + ad + be + bd + cd)q n 

-q n -\l-q n -q n+1 )abcd)), (9) 



9n{X) 



and y/5 = j-p 2 . The three-term recurrence f n+ i(x) = (x — a n )f n (x) — b n f n -i(x) for the monic 

Askey- Wilson polynomial can be found in the literature, for example eq.(3.1.5) in jTTj. Then 

^^Whilc preparing our manuscript, we became aware of a paper by van Diejen |18|. in which the same 
result was presented with a more elegant proof. 

2 This form of g(x) is obtained by using some empirical knowledge based on numerical analysis. 
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eq.(A.7) in [I] holds because the functions X n (x),Y n (x) of (A. 13), (A. 14) in [T] vanish for 

these a n ,b n and g n (x). Therefore we obtain Proposition A. 3, A. 4 in [Tj. 

Remark 1: The models studied in pQ are special cases of this model. For example, 

5eq.(2.75) in [T] (g , g u g 2 , g 3 , g 4 ) = (g L , \g s , \g s , 0, 0), 
C eq.(2.78) in Ij (g Q , g u g 2 , g 3 , g 4 ) = (gs, \gu \gu \ql, \ql), 
B'C eq.(2.79) in lj (g , g u g 2 , g 3 , g 4 ) = (g M , gs, gL, gL, 0), 

due to the following trigonometric formulas, 

sm(x — ig) \ cos(x — ig) \ 

= cosh gil — % tanh (7 cot x), = cosh g ( 1 + % tanh g tana;), 

sinx cos a; 

, n /H . , „ v sin(2a; — Via) sin(x — iq) cos(x — iq) 

cosh 2^(1 - « tanh 2^ cot 2x) = — ^ — = — — — . 

sin 2x sin x cos x 

Remark 2: Ismail et.al. studied the g-Strum Liouville problems and the Bethe ansatz equa- 
tion of the XXZ model ^E]- A special case of their results states that the zeros of Askey- 
Wilson polynomial p n (x; a, b, c, d\q) oc nj=i( x — cos ^Qj) whh the parameters 

q = e - 290 , (a,b,c,d) = ( e - 29l ,e- 292 ,e- 293 , e - 294 ), (10) 

satisfies the Bethe ansatz like equation © with (gj > 0) 

sin(x - ig ) sin(x - ig x ) sin(x - ig 2 ) sin(x - ig 3 ) sm(x - ig 4 ) 

v(x) — , w(x) = . (11) 

sin a; sinx sin a; sin a; sinx 

However this does not mean that {qj} are equilibrium positions of the system with (jllj) . 
because Vj(q) is not positive in this case. Moreover the cot 4 a; term in w(x) (fTTjl appears too 
singular to give a satisfactory quantum Hamiltonian. 

3 Single Particle Systems: Shape Invariance 

Next let us consider the shape invariance of the "discrete" quantum mechanical single particle 
systems with a g-shift type kinetic term. The argument is parallel to that given in in 
which discrete quantum systems with a shift-type kinetic term are discussed. 
In this section we use variables 6, x and z, which are related as 

0<6<n, x = cos6, z = e ie . (12) 
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The dynamical variable is 29 and the inner product is (f(9),g(9)) = d9f(9)*g(9). We 
denote D = D z d = z£. Then q D is a g-shift operator, q D f(z) = f(qz). Note that 

For a real constant g (0 < q < 1) and a function V(z) = V (z; A, g) with a set of real 
parameters A, let us consider the following Hamiltonian H = H(z; A, q), 



H^^y^)q D VV^+^VV^q- D VW)-Uy^+V*(^ ^) 
The eigenvalue equation reads 

H(f) n = S n <f> n , (15) 

with eigenf unctions <f> n (z) = (j) n (z;\,q) and eigenvalues £ n = S n (\,q) (n = 0,1,...) (we 
assume non-degeneracy So < £± < ■ ■ ■ ). The kinetic term causes a g-shift in the variable z. 
This Hamiltonian is factorized and consequently positive semi-definite: 

H = A*A, (16) 

A = A(z; A, q) d = f ( ? t v^^) - y/V{z)) , (17) 

^ = A(z; A,g)t ^ -L^v^i)^ _ v^pjp), (18) 

where f denotes the hermitian conjugation with respect to the above inner product. The 
ground state 0o is the function annihilated by A: 

A/> = (^#0 O = O, £ = 0). (19) 

Explicitly this equation reads 



sjv^q-^z^^z) = \Jv(q-h) o (<rM- (20) 
The other eigenfunctions can be obtained in the form 

<t> n {z) ^P n (z)Mz), (21) 

where P n ( z ) = P^ q \ z ) is a Laurent polynomial in z (for the explicit example below, it is 
the Askey- Wilson polynomial in x). This P n (z) satisfies 

HP n = S n P n . (22) 



Here H = H(z; A, q) is a similarity transformed Hamiltonian in terms of the ground state 
wavef unction </> : 

H ^ 4> 1 o H o O = §y(*) ^ + iy^^ 1 ) g- D - \{V{z) + (23) 

Corresponding to the factorization of H (jl6j) . is also factorized: 

H = BC, (24) 
C7 = C(«,g) = |(gT_ g -f), (25) 
B = B(z; A, q) = V(z) q% - V*(z~ l ) q~% . (26) 

Let us define a new set of wavefunctions <fii, n (z) = <f>i,n(z] A, q), 

0i, n = A/> n (n=l,2,...). (27) 

As a consequence of the factorization, they form eigenfunctions of a new Hamiltonian H\ = 
Hi(z;X,q), 

#i = A4 f (28) 

with the same eigenvalues {£„}: 

= AtfA<l> n = AS n( p n = £ n <j> ltn (n = 1, 2, . . .). (29) 

To consider the shape invariance of H, we try to find the operators A 1 , A\ and a real constant 
£i satisfying 

H x — AA^ — AXA-i + E x , (30) 
A x = A x {z\ A, g) d ^ -±= (qt VWi^ 1 ) ~ <T f \ZW)) , (31) 

A\ = A 1 (x; A, g)t d ^ f -L^^V^qf - y/V^j ) . (32) 
In other words, given V(z) = V(z; A, q), find a new potential Vi(z) = V\(z; A, q) satisfying 

Vx^V^q^z- 1 ) = V(qh)V*(q-h- 1 ), (33) 
V 1 {z) + V*(z~ l ) = V(q-h) + V*(q-h- 1 ) + 2Sl (34) 

If Vi has the same functional form as V with another set of parameters A' (e.g, g-shifted A), 

Vi(z; A, q) oc V(z; A', q), (35) 



then it is shape invariant. Suppose V\ has the form 

V 1 (z) = V(qh)g(z), (36) 

with an as yet unspecified function g(z), the above conditions (jSHJ), ffi^jl get slightly simpli- 
fied: 

g(z)g*(q- 1 z- 1 ) = 1, (37) 
Vfas 3)0(2) + V^z -1 )^* -1 ) = V(q-h) + V*(q-h~ 1 ) + 2£ v (38) 

If the desired V\ is found, we can construct H 2 ,H^,--- by repeating the same step. We 
illustrate this procedure by taking the Askey- Wilson polynomial as an example. 
Let us take V as 

"M - n*, a, «) - (1 "'^:^ 1 :g ) (1 " fe) . (^) 

where A = (a, 6, c, <f). For simplicity we assume — 1 < a, b, c, d < 1. Note that V*(2) = V(^). 
The ground state (fTHj) is given by [TT] 



<f>o(z) = 4> {z- } X, q) ocy/w(z; X, q) 



def 



(z 2 ;q) 



00 



(az, bz, cz, dz; q) c 



(40) 



y z 2 , z 2 , 5)00 (41) 



(az,az 1 ,bz,bz l ,cz,cz 1 ,dz,dz 1 ;q) c 



where ( a\ , a m 

j 9)00 = njli Il^Lo(^ — a j < f 1 )- Excited states have the form (|2*T|) n (z) oc 
-f'n(^)0o(- 2 ), and implies that P n {z) is proportional to the Askey- Wilson polynomial [TT], 

P n (» = Fi A ' 9) (^) oc p n (x; a, b, c, d; q), (42) 

E n = £ n (X, q) = \q- n (l - q n ){\ - abcdq^ 1 ), (43) 

which is an orthogonal polynomial, 

w(z; X, q) p n (x; a, b, c, d\q) p m (x; a, b, c, d\q) oc S nm , (44) 



1 VT^ 2 

namely (4> n ,(j>m) oc 5 nm . By denoting p n (x; a, b, c, d\q) = P n (z;X,q), the factorization (|2*4l) 
gives the forward and backward shift relations ((3.1.8) and (3.1.10) in [TTj): 

(7(2;; q)P n (z; X, q) = -S n (X,q)q^(z - 2 _1 )P n _i(z; q^A, g), (45) 
-B(z;X,q) q^(z- z- l )P n _ 1 (z;qiX,q) = P n (z;X,q). (46) 
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It is easy to check that V\ in the form (J35)l with 



1 - a 2 z 2 

9{z) = <f X \^T ( 4 7) 



satisfies (}3Tj) and (J35j) . and it becomes 

V x {z- X, q) = V(qh)g(z) = gU, g), (48) 

and £i is 

5 1 (A,g) = ^g- 1 (l-g)(l-o6cd). (49) 



Therefore we have shape invariance, (J48|) and 

A 1 {z;\,q) = q-*A{z',q*\ } q), (50) 

H^z; X, q) = q-'Hiz; gU, q) + £i(A, g). (51) 

We write down important formulas once again: 

A(z;X,q)M^\q) = 0, (52) 

A{z- A, g)^; A, g)t = g" 1 ^; g*A, g^A, g) + ^(A, g). (53) 

Starting from V = V , H = H , Oi „ = <p n , let us define V s , H s , (p s ^ n (n > s > 0) step by 
step: 

V s+1 (z;X,q)^q- 1 V s (z;q 1 iX,q), (54) 

# s+ i(z; A, q) d = A s (z; A, g)A s (^; A, g) 1 " + £ S (A, g), (55) 

(p s+hn (z; A, g) = f A s (z; A, q)<j) s ^{z] A, g). (56) 

Here yl s and are defined by 

A s (z; A, g) d = f -L ( g # ^V/^-i; A,g) - <T VVto A,?)) , (57) 

A s (z; A,g)t ^ -L(VK(0;A,g)gt - ^/(z" 1 ; A, g) g~#) . (58) 
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As a consequence of the shape invariance (J33J), we obtain for n > s > 0, 

V s (z;\,q) = q- s V(z;qi\,q), (59) 

A s (z; A,g) = q-iA(z;q*\,q), A s (z;\q) ] = q~i A(z;qi\q)\ (60) 
H s {z- A, q) = A s {z; A, qfA„{z\ A, g) + £ S (A, q) = q~ s H{z; q*\, q) + £ S (X, q), (61) 

£ s+ i(A,g) =£ s {\,q)+q- s £ 1 {qi\,q), (62) 

# s (x; A, g)0 s ,„(2;; A, g) = £ n (X, q)(j) s , n (z; A, g), (63) 

4 a (z; A,g)0 s , s (z; A, g) = 0, (64) 

A s (x; A, g)Vs+i,n(^; A, g) = (£„(A, g) - £ S (A, q))(f> Stn (z] A, g). (65) 

The relation (Jfi2|) means that {£ n } is calculable from £\ (I49j) . In other words, the spectrum 
is determined by the shape invariance. 
From (J56j) and (J65|) we obtain formulas, 

(f) S:n (z; A, g) = A s ^(z; A, g) • ■ ■ A, q)A (z; A, g)</> n (2; A, g), (66) 

, / x ^ (z;A,g) t A,g) f A n _ x {z\ A, g) f 

gj " £ n (A, g) - £ (X, q) £ n (X, q) - £,(X, g) " ' £ n (X, g) - ^(A, g) A ' q) ' 

(67) 

The former gives the eigenfunction Sjn of the s-th Hamiltonian ff s along the isospectral 
line with energy £ n , starting from (j) n of the original Hamiltonian H by repeated application 
of the A operators. The latter ([67)1. on the other hand, expresses the n-th eigenfunction 
(j) n of the original Hamiltonian, starting from the explicitly known ground state 4>n,n of 
the n-th Hamiltonian H n by repeated application of the A^ operators. Since (j6ip implies 
4>n t n(z; A, q) oc 4>o(z; q^X, q), (j) n is expressed in terms of O an d V. The latter formula (|57jl 
could also be understood as the generic form of the Rodrigue's formula for the orthogonal 
polynomials. The situation is depicted in Fig.l. The operator A acts to the right and A^ 
to the left along the horizontal {isospectral) line. They should not be confused with the 
annihilation and creation operators, which act along the vertical line of a given Hamiltonian 
H s going from one energy level £ n to another £ n ±±. 

In order to define the annihilation and creation operators, let us introduce normalized ba- 
sis {(p s ,n}n>s for each Hamiltonian H s . Ordinarily, the phase of each element of an orthonor- 
mal basis could be completely arbitrary. In the present case, however, the eigenf unctions 
are orthogonal polynomials. That is, they are real and the relations among different degree 
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£i 03 
( 

£-2 4>2 

( 

£1 0i 
£ 

00 

Ho Hi H-i H 3 

X qi\ qi\ qi\ 

Figure 1: A schematic diagram of the energy levels and the associated Hamiltonian sys- 
tems together with the definition of the A and A^ operators and the 'creation' (aJ) and 
'annihilation' (a) operators. The parameter set is indicated below each Hamiltonian. 

members are governed by the three-term recurrence relations. So the phases of {<p s ,n}n>s 
are fixed. Let us introduce a unitary (in fact an orthogonal) operator U s mapping the s- 
th orthonormal basis {4> s ,n}n>s to the (s + l)-th {4> s +i,n}n>s+i (see Fig. 1 and for example 

nana): 

U s <Ps,n = 4> s+ l,n+l, U^4>s+l, n +l = 4>s,n- (68) 

We denote that U = U. Roughly speaking U changes the parameters from A to q^X. Let 
us introduce an annihilation a and a creation operator for the Hamiltonian H as follows: 

a = a(z; A, q) d = U^A(z; A, q), a f = d(z; A, g) 1 " = f A(z; A, q)^U. (69) 

It is straightforward to derive 

H = a f fi, (70) 
[a, a)]4> n (z\ A, q) = (£ n+1 (\, q) - £ n (X, q))(j) n (z] A, q). (71) 

4 Summary and Comments 

In this article we have studied the equilibrium positions of the Ruijsenaars-Schneider-van 
Diejen systems with the trigonometric potential and shown that for a suitable choice of the 
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elementary potential functions they are given by the zeros of the Askey- Wilson polynomials 
with five parameters (see the footnote in section^. The equation for the equilibrium posi- 
tions (JHJ) (without the positivity condition) can be written in the Bethe ansatz like equation 
(jHJ. The Bethe ansatz is a powerful method for solvable models, and solving the Bethe 
ansatz equation or clarifying the properties of its solutions are very important. Ismail et.al. 
studied Bethe ansatz equations for spin s XXZ models from the g-Sturm-Liouville problem 
point of view [TSj. This kind of approach would shed new light on the Bethe ansatz(-like) 
equations. 

We have also studied the shape invariance of "discrete" quantum mechanical single par- 
ticle systems with a g-shift type kinetic term. As an example of this shape invariance, we 
present such a system whose eigenfunctions are the Askey- Wilson polynomials. In this ex- 
ample Viz) is a rational function of z (trigonometric function of 6), but the method works 
for a wider class of functions. In ordinary quantum mechanics there is the Crum's theorem 
[T3] . which states a construction of the associated isospectral Hamiltonians H s and their 
eigenfunctions (Fig.l) even if the system has no shape invariance. The construction of 
H s and (p S)Tl given in this article and 3] needs shape invariance. A "discrete" analogue of 
the Crum's theorem, namely similar construction without shape invariance, would be very 
helpful, if exists. 

We comment on the shape invariance of the Askey- Wilson polynomials with a small 
number of parameters, p„(x; a, b, 1, —1). In this case V(z) in (jT4*)l is 

„, , N (l — az)(l — bz) . , 7N 

V{z\\q) = - r^- 2 S A=(a,6). (72) 

1 — qz z 

By taking a form fpTFTf> with the same g{z) (jUJ), the conditions (J37)) and (jHEJ) are satisfied, 
and we have 

V 1 {z;\q) = q- 1 V(q-h;q\,q), £ x i\q) = ^(q' 1 - l)(l + o6). (73) 

Since D z = z-^ is invariant under the rescaling of z, D z = D az , the Hamiltonian is shape 
invariant, 

#i(3,A,g) =q- 1 H(q-h;q\,q)+£ 1 (\,q). (74) 
The s-th Hamiltonian and the spectrum are given by 

H s {z; A, q) = q- s H(q-i Z] q s X, q) + E s {\, q), (75) 
S a {\ q) = £ S _!(A, q) + q-^e^X, q) = l -q~ s {l - q s )(l + abq^ 1 ). (76) 
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